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Abstract The frequency-moment expansion method is developed to analyze the validity of the Luttinger sum rule 
within the Mott-Hubbard insulator, as represented by the generalized Hubbard model at half filling and large U. For 
the particular case of the Hubbard model with nearest-neighbor hopping on a triangular lattice lacking the particle-hole 
symmetry results reveal substantial violation of the sum rule. 

PACS. 71.10.-w Theories and models of many-electron systems - 71.27.-l-a Strongly correlated electron systems; 
heavy fermions - 71.18.-l-y Fermi surface: calculations and measurements; effective mass, g factor 



1 Introduction 

The concept of a Fermi liquid (FL) is the basic one for the un- 
derstanding of electrons in solid state. The central pillar of the 
FL theory is the existence of a Fermi surface with a volume 
unchanged by electron-electron interactions, having the firm 
support in the Luttinger sum rule (LSR) |1,2, 3|. In last two 
decades experiments on novel materials with strongly corre- 
lated electrons, in particular superconducting cuprates [4 |, and 
theoretical analyzes of corresponding microscopic electronic 
models revived the question of the possible limitations and the 
breakdown of the LSR. 

The LSR within a paramagnetic metal relates the density 
of electrons to the Fermi surface (k space) volume defined by 
poles of the Green's function (GF) G(k, lo = 0). Only recently 
it has been recognized that the concept can be generalized to 
insulators |5,6| whereby the coiTesponding 'Luttinger' surface 
is defined by the zeros of G(k, 0) — 0. It has been shown that 
the LSR is indeed satisfied for models with a particle-hole sym- 
metry at half filling, and in particular for the frequently invoked 
Hubbard model on a one-dimensional (ID) chain and on a 2D 
square lattice |7|. There are several indications that in the ab- 
sence of particle-hole symmetry the LSR can be generally vio- 
lated |l7][8_| although no explicit example emerging from a mi- 
croscopic model with a repulsive Coulomb interaction term and 
within the canonical ensemble (assuming fixed particle den- 
sity) has been presented so far. On the other hand, such cases 
seem to be found for a Hubbard model on an inhomogeneous 
lattice and on a homogeneous lattice with odd number of sites 
close to half filling |9|. Recently, the present authors pointed 
out on the generalization of the LSR to finite systems of inter- 
acting electrons I lOj . Since the original proof [2 1] of the LSR 
remain valid even for small systems, this allows an alternative 
way to show and understand a possible breakdown of LSR. So 



far, results confirm cases of violation in restricted (nonpertur- 
bative) models as the t-J model, whereas reachable systems 
and discussed model regimes did not provide a clearcut viola- 
tion within the translationally symmetric Hubbard model | lOl. 

In this paper we analyze in more detail the validity of the 
LSR within the insulating state emerging from a single-band 
metal via the Mott-Hubbard mechanism of strong electron-ele- 
ctron repulsion. As the prototype model we use the Hubbard 
model on triangular lattice with the large onsite repulsion U ^ 
t. We develop an approach employing the expansion of fre- 
quency moments in powers of t/U to determine the value of 
the G{k, 0) leading to the location of the Luttinger surface 
G(ki, 0) — 0. The method is valid in the thermodynamic limit 
(for an infinite system) but can be tested as well on small sys- 
tems, which can provide also required static quantities as in- 
put. While the method confirms the validity of LSR in Hub- 
bard models with a particle-hole symmetry [7 1, we further con- 
centrate on the Hubbard model on a triangular lattice as the 
case without the latter symmetry. Results show that the LSR 
is violated substantially for large U /t, being an indication for 
analogous violation in more general Mott-Hubbard insulators 
without particle-hole symmetry. 

For simplicity we restrict in the following our discussion to 
the single-band Hubbard model on a general lattice 



H 



ijs 



(1) 



where U is the local Coulomb repulsion, while as the inter- 
site hopping tij we can consider besides the nearest neighbor 
(n.n.) term tij = t also extended model, e.g., with the next- 
nearest neighbor (n.n.n.) hopping tij = t' . We are interested 
in the regime of large enough U ^ tij and half filling n = 1, 
where the ground state is insulating with a Mott-Hubbard gap 
separating both Hubbard bands. Reduced models, such as the 
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t-J model, are mostly adequate enough to discuss physics in 
this situation. The essential motivation to stay with the Hub- 
bard model is the fact that the latter model is perturbative, at 
least for small U, and thus could be adiabatically connected to 
the case of noninteracting (U = 0) fermions. The breakdown 
of the LSR could be then attributed to nonadiabatic develop- 
ment of the ground state (i.e. in the limit T ^ 0) through the 
metal-insulator transition, allthough there could be also other 
interpretations. 

In the following we first derive lowest orders of GF with 
the use of large U expansion method. We also use particle- 
hole symmetry and give in some more detail the calculation 
of chemical potential. In chapter[3]we first give general results 
of large U expansion method and in chapter 13.11 we present 
our results for particular case of triangular lattice. At the end a 
discussion of results and conclusions are given. 



Analogous procedure can be repeated for G+ with the result 



G+(k,0) = ^(- 



ri=0 



2 s. n+1 

Tj' 



M+(k) = (ck.[i?',...[i/',4j]), (7) 

where H' = H — U{Ne ~ N), with being electron-number 
operator and N is the number of sites in the system. G+ is also 
tightly related to with the particle-hole transformation, as 
discussed lateron. 

Let us now present first two moments explicitly. 



-ik(i-j) / t 



N 



2 Large U expansion 

The central quantity for the LSR is the GF at temperature T 
0, given by 

G{k,uj) = -i r dte^('^+^)*(0|cLck.W|0) 



+1 / dte'(-+^)*(0lcks(t)4j0), (2) 



where /i is the chemical potential. Representing the GF in terms 
of the spectral function 



G(k,w) 



-duj'. 



(3) 



We deal with the insulator with real G(k, oj ^ 0) and a large 
gap A ^ U separating regions of nonzero spectral functions 
yl(k, ut). In order to evaluate the relevant G(k, 0) we separate 
contributions 



G(k,0) = G-(k,0) + G+(k,0) = 



-oo M - ^ - 



^ — uj' + if] 



du' (4) 



Since /i ~ C//2 it is reasonable to expand both G (k, 0) and 
G+(k, 0) in powers of {Uj/U)". 
Let us first consider G^ , 



. n "-^ c 



A(k,Cj'){Cj' - fL)"d(:;' 



n=0 

oo 



= E(|)""'E^^^^™(k)r)(-Ar, (5) 

ri=0 m=0 ^ ^ 

where jl ~ fjL — U /2. Note that Mj^ (k) are frequency moments 
which can be directly expressed in terms of the equal-time ex- 
pectation values within the half-filled ground state |0) involv- 
ing I commutators with H, 



Mr(k)-([ff,...[i/,4j]cks). 



(6) 



N 



where e(k) is the free band (U — 0) dispersion as determined 
by hopping parameters tij . 

One expects that for large U ^ ty , is of the order of 
t\j, at least for lowest orders I, as well as /i cx tij . Hence the se- 
ries Eq.Q should be well converging. To show that Mj^ oc t\p 
we employ at this stage the well known canonical transforma- 
tion of the Hubbard model for large U [11,12 ,131, with which 
the number of doubly occupied sites becomes a good quantum 
number. This leads at half-filling, n = 1, to an effective spin 
model fT2l and at finite doping to strong coupling (SC) model 
for the lowest orders ||T3]| . 



6 c 6 



(9) 



Within the lowest order of the tij /U expansion one gets S = 
= -jjYl K5cLcjs(l-nj5) -il-n,s)clcjsnjs]. (10) 

ijs 

Transforming expressions for moments (O and using Eq. 
([Tol l, we arrive to explicit expressions at n = 1, 

M^ik) = 1 ±|5]e,(k)KS, . So)-i] + 0(^), 



U 



Mf (k) = ^ + ^£,(k)[(S, . So)-i] + 0(^), (11) 



where es{k) refers to partial 'free' bands corresponding to par- 
ticular neighbors (n.n., n.n.n. etc.), respectively 



£5(k) = -i5^e'''---^ 



(12) 



We note that spin correlations {Ss • So) in Eq. (fTTT l are to be 
evaluated in the ground state of iJg/ and that such state does 
not contain doubly occupied sites. 
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2.1 Particle-hole symmetry 

Let us now return to G+. In order to show that its moments 
are closely related to Mj~, it is very instructive to realize that 
can be obtained from G~ with particle-hole transforma- 
tion, which changes Cis c]^ and vice versa. Applying the 
transformation to in Eq. Q we note that also the ground 
state |0) changes to |0), corresponding to transformed H. H is 
obtained from H by replacing tij — > —tij and adding the extra 
potential term U {N — Ne). We therefore arrive to the following 
expression for G+, 

G+(k,0)=^/" dte*(^"'')*(0|4^e'^*Ck.e-*-f^*|0). (13) 

which is analogous to G^ in Eq.©. One can then continue by 
analogy to G~ using the unitary transformation with S, which 
is the same form as S with —tij instead of Uj . State, for which 
expectation value must be calculated, is now the ground state of 
Hef. Moreover, for half filling at n = 1, H^f contains only spin 
operators, whereby odd spin terms (also odd terms of ty) are 
forbidden by the isotropy of original Hamiltonian, and /fe/ = 
Hef lfT2l . and therefore also ground states of them are the same. 
Finally we can relate 

G+(k,0)--G-(k,0)^^_^, (14) 

where G^ stands for G^ with tij changed to —tij and /i 
—fi. Eq. (fT4l) confirms that Mj^ are simply related to Alf' as 
already seen from relations (fill) . 

Expressions are so far valid for general hopping tij, with 
the only assumptions that we are dealing with the paramagnetic 
state with total spin 5 = and without broken translational 
symmetry. 

2.2 Chemical potential 

We are interested in the case of Mott-Hubbard insulator at half 
filling, n=l and within the canonical ensemble the chemical po- 
tential /i should be determined from the grand partition func- 
tion by first fixing the number of particles in the system and 
then taking the limit T — > 0. Such a procedure gives an unique 
value for /i in the middle of the gap |8,14,15|, and works as 
well for finite systems 1 10|. Different limits with /i being any- 
where within the Mott-Hubbard gap have been also recently 
considered |l8][9l. Our procedure simply locates fi (as well as in 
any finite system with A^e particles 1 10|) with the ground state 
energies of systems with one electron more and less, respec- 
tively iiiMKia, 

M=(i?o^=+^-<-^)/2 (15) 

where for half-filling = N. Since the ground state of one 
hole for Hf^f does not contain doubly occupied sites, we get 
Eq~^ = 0{tij). On the other hand, ground state with one 
added electron more contains one doubly occupied site and 
therefore Eq^^ = U + 0{tij). fj, is therefore within the gap 
and approximately /i ~ C//2. Energies E^^^, E^^^ as well as 
ji might be expanded analytically in terms of tij/U, however 
we use furtheron numerical results instead. 



3 Results 

Using expanded moments (fTTI) in G^ ^ and Eq. (fT4l l for G+, 
we finally get for the total GF Q, 

G(k,0) = 4£5(k)(S, . So) - A) + O(^), (16) 

s 

From Eq. ^ and (fT4l i one can also conclude, that for fl = 0, 
GF consists only of odd terms in tij . 

It has been already realized Q, that the LSR remains valid 
for cases with particle-hole symmetry, where H = H. This 
follows also from our analysis. If we assume in Eq. ( fTSI l only 
nearest neighbors and put p, — Q, which is the consequence 
of the particle-hole symmetry, GF changes sign at same k as 
e(k). Within the lowest order in tij/U'^, this leads to the same 
Luttinger volume as for noninteracting fermions. 

In order to find cases where the LSR might be violated, we 
turn to systems without the particle-hole symmetry. Here, in 
general /i ^ and also £(k) is not negative for exactly half of 
the Brillouin zone. Then, Eq. ( fTSI l seem to suggest that the LSR 
must be violated in most cases, since the perfect cancellation of 
two quite different terms in Eq. ( fT6] l looks quite improbable. In 
the following, we show that such compensation indeed does not 
occur for one of simplest non-symmetric cases, i.e. within the 
Hubbard model with only n.n. hopping on the triangular lattice. 

3.1 Hubbard model on a triangular lattice 

Furtheron we consider the model with n.n. hopping only, tij = 
t. To check the LSR on triangular lattice at least to lowest order 
in t/U, some quantities are needed. Free band dispersion e(k) 
is given with Eq. ( fT2] l and equals ei(k). Spin correlations Cs = 
{Ss ■ So) should be evaluated in the ground state of H^f at 
half-filling. Within the first order in t'^ /U this corresponds just 
to the Heisenberg model on a triangular lattice. The value for 
n.n. spin correlations in Heisenberg model has been evaluated 
by a number of authors lT6l[T7l[T81 and Gi = -0.182. This 
value is quite insensitive to finite t'^/U, e.g. we obtain within 
the Hubbard model on = 15 sites -Gi - 0.183 - 0.189 for 
tVC/ = 0.05 i - 0.02 t. 

The chemical potential fl is more delicate quantity. Here, 
we estimate jl as accurate as possible by calculating Eq^^ 
and E^^^, performing the exact diagonalization of finite clus- 
ters using the Lanczos algorithm. Lattices with periodic bound- 
ary conditions are chosen in such way, that they allow three- 
sublattice order 1 16] and have a minimal imperfection 1 19 1, that 
is zero for all choosen lattices, except for N=18 it is 2. The 
maximum size for the Hubbard model is = 15, for which 
we present /2 in Fig.[T] Fortunately, for larger U jt the Hubbard 
model can be well represented by the SC model (equivalnet to 
the t-J model additionally including three-site correlated hop- 
ping term of the order t^ /U) llT3l . As seen in Fig. [T] jl cal- 
culated from the Hubbard model and from the SC model for 
TV = 15 match well for small t^ /U < 0.05. This gives support 
to the values obtained for fl within the SC model, as presented 
in Fig.[T]for systems up to iV = 24. 

In considering fl one should take into account that at large 
U/t the Nagaoka ferromagnetic state is possible, i.e., E^^^ 
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Figure 1. (Color online) Chemical potential /i on a triangular lattice 
vs. t^/f/ calculated on clusters with TV = 15 — 24 sites within the 
Hubbard model (15Hub) and in strong coupling (SC) models, respec- 
tively. The value jli for which LSR would be satisfied is also shown. 



Table 1. Numerically obtained expansion coefficients. 



k 


B(k) 


(0,0) 


-2320 ±40 


(27r/3,0) 


-167 ±3 




200 ±5 


(47r/3,0) 


1810 ± 40 



corresponds to 5* = N/2 or at least not to 5 = 0, 1/2. In Fig.[T] 
we therefore present results only within the regimes of lowest 
S. In finite systems the onset of Nagaoka state shows a size 
dependence and moves to lower values P/U with increasing 
sizes. From our analysis we can estimate, that the paramagnetic 
state with S = (for even + 1) remains stable at least for 
P/U > 0.025i. Taking mentioned limitations into account, 
we can analyse results in Fig. [T] which give p, ^ — 0.19t + 
6.8t'^/U ± O.lt for 0.025t < t'^/U < 0.051 Unfortunately, 
the proper finite size scaling cannot be performed since lattice 
shapes as well as S are not the same for different N therefore 
the uncertainty is an estimate. 

Before presenting the final result, we comment on the higher 
order corrections to Eq. ( fT6b on a triangular lattice. With the 
use of particle-hole transformation we can show that for fi ~ 
there is no correction of order /U^ and the next correction of 
order t^/lf^ can be presented as 



G(k,0)^. 



16 

U 



So) + B{k) 



(17) 



We evaluate -B(k) for some k with the use of numerical cal- 
culations on = 12 site Hubbard cluster -B(k) is obtained 
by the expansion of numerically calculated GF for the range 
U = 60t — 120i. Results are presented in Table [T] A relevant 
observation is that the higher order correction is small in vicin- 
ity of the Luttinger surface for large enough U > 20 t. More- 
over, the correction is negative and suggests (even) smaller Lut- 
tinger volume, which is the tendency towards larger LSR vio- 
lation. 




-0.006 



0.012 



Figure 2. (Color online) Green function G(k, 0) for U/t = 40 within 
the first Brilloiun zone. Displayed are also the Luttinger volume, 
which would satisfy Luttinger theorem (dashed line) and calculated 
Luttinger volume (full line). 



Our final result for GF can be thus presented as 



G(k,0) 



4 

IP 



(4Cie(k) - G.19i + 6.8f/U ± O.lt) . (18) 



In Fig.|2]we show results, Eq. ( fTSl l. obtained with the expan- 
sion method for particular U — 40t. Dashed black line shows 
the Luttinger volume, which would correspond to noninteract- 
ing fermions. The latter can be obtained from Eq. ( fTSl ) if one 
would insert pi = —O.Qlt instead of the calculated value. Evi- 
dently, pi is far from any reasonable estimate in Fig.[Tl where 
fi{U = 4Gi) 0. Therefore the calculated Luttinger volume 
is obviously smaller than required by the LSR, the difference 
being about 17% in Fig.|2l Estimated uncertainty is too small 
to compensate the violation of the LSR. 



4 Conclusions 

In conclusion, we presented evidence that the LSR is violated 
on triangular lattice for half-filled Hubbard band in regime of 
large U /t. Clearly we should stress again limitations to our 
analysis and results. First, U /t must be large enough so that the 
ground state corresponds to an insulator with a Mott-Hubbard 
gap which is the case for the Hubbard model on triangular lat- 
tice at U ^ 12t |20|. On the other side, at very large U/t the 
onset of Nagaoka instability appears for the 7V+ 1 ground state. 
In absense of other studies we estimate that this can happen 
only for U > 40t. 

Finally, let us speculate on the origin of the violation of 
the LSR. A frequent mechanism for the latter is related to the 
breakdown of the paramagnetic state or the loss of translational 
symmetry, in a microscopic system connected with a phase 
transition and the onset of the coiTesponding long-range or- 
der Besides discussed Nagaoka instabilities possible at very 
large U/t, within the triangular lattice another possibility is the 
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long-range antiferromagnetic order which has been established 
within the Heisenberg model 1 16] but so far not within the Hub- 
bard model at half filling |21 1. The onset of such order would 
limit the relevant U /t parameter window of the paramagnetic 
insulating state. On the other hand, any breaking of transla- 
tional symmetry is expected to change entirely the Fermi sur- 
face topology and not to partially reduce its volume as found in 
our study. 

More radical conclusion would be that the LSR is gener- 
ally violated within the Mott-Hubbard insulators without the 
particle-hole symmetry |7 ,8|, and in particular within the in- 
sulating state of the Hubbard model on a triangular lattice. We 
note that the main formal requirement for the validity of the 
LSR im is the adiabatic development of the ground state (i.e. of 
the free energy in the limit T ^ 0, where the LSR becomes ap- 
plicable) with the increasing electron-electron repulsion. Hub- 
bard model is perturbative in U and adiabatically connected to 
the noninteracting fermion system at least for weak U /t. On 
the other hand, the metal-insulator transition at t/ — Uc seems 
to represent a point, where the perturbation theory as well as 
LSR breaks down. Such scenario is still far from obvious so 
further studies in this direction are needed. 
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